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We investigate the doping of AA-stacked graphene bilayers. Applying a mean field theory at zero 
temperature we find that, at half-filling, the bilayer is an antiferromagnetic insulator. Upon doping, 
the homogeneous phase becomes unstable with respect to phase separation. The separated phases 
are undoped antiferromagnetic insulator and metal with a non-zero concentration of charge carriers. 
At sufficiently high doping, the insulating areas shrink and disappear, and the system becomes a 
homogeneous metal. The conductivity changes drastically upon doping, so the bilayer may be used 
as a switch in electronic devices. The effects of finite temperature are also discussed. 
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Introduction. — Controlled metal-insulator (M-I) tran- 
sitions are a very useful property for electronic appli- 
cations of graphene Such transitions have been ana- 
lyzed theoretically (e.g., and experimentally observed 
in graphene by several groups using different techniques, 
e.g., chemical adsorption H, thermal annealing gate- 
induced M-I transition [5|, and percolation-driven M-I 
transition in graphene nanoribbons due to inhomoge- 
neous electron- hole puddle formation . 

Here we study an AA-stacked bilayer of graphene ( AA- 
BLG). The purpose of this work is to demonstrate that 
this systern, which has been recently successfully fab- 
ricated 0, Q, can exhibit a M-I transition upon dop- 
ing. Further, we will demonstrate that the required lev- 
els of doping are within current experimental capabili- 
ties. Unlike AB-stacked bilayers, the AA-BLG received 
very modest theoretical attention However, ad- 

vances in fabrication of AA-stacked bilayers and multi- 
layers 0, [1] underscore the need for thorough theoretical 
investigations. 

Tight-binding calculations for AA-BLG d [l^ predict 
that near the Fermi energy the bilayer has two bands, one 
electron-like and one hole-like. These bands have Fermi 
surfaces, unlike Fermi points in monolayer graphene and 
AB-stacked bilayers. An important feature of the AA- 
BLG is that the hole and electron Fermi surfaces coin- 
cide. As shown in Ref. Il3l . if interactions are included, 
these degenerate Fermi surfaces become unstable, and 
the bilayer turns into an antiferromagnetic (AFM) in- 
sulator with a finite gap. This electronic instability is 
strongest when the bands cross at the Fermi energy. Im- 
purities or doping shift the Fermi level and suppress the 
AFM instability. 

Superficially, one may expect that the AFM gap A de- 
creases with doping x and vanishes above some critical 
value Xc- However, we will show that the homogeneously- 
doped state is unstable with respect to the phase separa- 



tion into undoped AFM insulator and doped metal. As 
the doping grows, the concentration of the AFM insula- 
tor shrinks, while it grows for the metal. Above a cer- 
tain threshold x* , metallic islands connect into an infinite 
cluster, and the percolation-driven insulator-metal tran- 
sition occurs, at which point the sample becomes metal- 
lic. 

Here we study the electronic properties of the doped 
AA-BLG in the framework of the Hubbard-like model 
used in Ref. [l^- We determine how the gap A depends 
on X in the homogeneous state and find the critical value 
Xc, where A vanishes. We further show that at small 
doping the homogeneous state is unstable because the 
compressibility of the system is negative, and find the 
doping range where this instability arises. The effects of 
non-zero temperature are also discussed. 

The model. — The Hamiltonian for pz electrons of car- 
bon atoms for the AA-BLG can be written as 



H — Ho + Hir, 



(1) 



where Hq describes electron hopping and Hi^t is the 
electron-electron interaction, fi is the chemical potential, 
and N is the operator of the total electron number in the 
system. In the tight-binging approximation 



Hn 



— t y ^ ^nia^mi 
{nm)2(T 

\ Tier 



(2) 



H.c. 



Here aj^j^. and a^,^^ (^Licr ^"^^ ^mia) ''^^'^ ^^'^ creation and 
annihilation operators of an electron with spin a in the 
layer i = 1, 2 on the sublattice A {B) at site n £ A 
(m e B). The amplitude t {to) in Eq. ([2]) describes 
the in-plane (inter-plane) nearest-neighbor hopping. For 
calculations we will use the values of the hopping inte- 
grals t « 2.57eV, to ~ 0.36 eV specific to multilayer AA 
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FIG. 1: (Color online) The band structure for the homo- 
geneous phase of the AA-stacked bilayer graphene near the 
K point; k = K + Skyey. The dashed lines show the non- 
interacting single-electron bands. At half filling these bands 
intersect with each other at the Fermi energy fj, = 0. Adding 
interactions opens a gap. The mean field bands [see Eqs. (O] 
are shown by solid lines. With doping, these bands are filled 
up to the level n' — fi — Ux/2. As a result of doping, the Fermi 
surface degeneracy disappears, and we have two Fermi surface 
components around each Dirac point. The inset shows the 
first Brillouin zone (hexagon) and the reciprocal lattice unit 
cell (rhombus) of the AA-BLG. Circles around the K and K' 
points correspond to Fermi surfaces of the doped system. 



systems [1^. Longer-range hoppings are neglected be- 
cause these are small (about or less than 0.1 eV), and 
we checked that the effects they produce are negligible 
(within 1-2%). 

The on-site Coulomb interaction can be written as 



Hint — X/ y^ni-Aa - ^ 



1 



(3) 
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where n^Aa = O-maO-ma- IT-miBa 

—a. It is known that the on-site Coulomb interaction 
in graphene and other carbon systems is rather strong, 
but the estimates available in the literature vary consid- 
erably 15, 3], ranging from U ^ t to ^ At. Because of 
this uncertainty, we will present our results in the form 
of [/-dependent functions, rather than definite estimates. 

Antiferromagnetic state. — In the absence of electron- 
electron coupling, C/ = 0, and zero doping (a; = 0, which 
corresponds to half filling) the AA-BLG band structure 
is shown in Fig. [1] by dashed lines. Two bands pass 
through the Fermi energy level near the Dirac points 
K = 27^{^/3, l}/(3^/3a) and K' = 27r{v^, -l}/(3V3a), 



where a is the in-plane carbon-carbon distance. The 
chemical potential is /i = 0, while the Fermi surfaces are 
given by the equation |/k| = to/t, where 



/k = 1 + 2exp(3ifca;a/2)cos(fcyaV3/2 



(4) 



For to/t <^ 1, one can expand the function |/k| near the 
Dirac points and demonstrate that the Fermi surface con- 
sists of two circles with radius fc,. = 2tQ/{3ta) around the 
Dirac cones K and K'. Upon doping, these Fermi sur- 
faces arc transformed into four circles [sec the inset in 
Fig. [T]. The presence of two bands with identical Fermi 
surfaces makes the system unstable with respect to spon- 
taneous symmetry-breaking. 

Since the unit cell of AA-BLG consists of four 
atoms, it is convenient to introduce the bi-spinors 

"V^L = (v4.AfT' V'kB^)' '^ith spinor components V^^^ = 

{alia^al^a) and Vkg^ = e'"^"- {blia^ ^l^a) ^ where 
— arg{/k}. The Hamiltonian Hq in this basis is 
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(5) 



In mean- field, the interaction operator Hi^t, Eq- ©j 
is replaced by a single-particle operator which breaks a 
certain symmetry of the system. As it was shown in 
Ref. [13 the ground state of our model is G-type AFM 
(that is, the spins on any two nearest-neighbor sites are 
antiparallel) , for which the spin- up and spin-down elec- 
tron densities are redistributed as rti^^ = n2Bt — ^2Ai = 
niBi = (l+a;+An)/2 and uiai = "2i3; = n2At = "let = 
(1 + a; — An)/2, while the total on-site electron density 
n = Uiaa + niaa = 1 -|- x is the same for any site. The 
mean-field interaction Hamiltonian for such phase is 



n. 
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Ux 



^ + A ^ (^k.AJ.'^-^^kyti - i>l,A-\^ ^'^^M (6) 
k 

"V'kSi'^^^kBi + V'ket'^^^ket , 



where is the Pauli matrix, and A ~ is the 

AFM gap, which should be found self-consistently. 

To find the gap, we solve the corresponding Shrodinger 
equation and derive the expressions for four electron 
bands E"^ (k) and eigenvectors 



(7) 



(2,3) 



TVA2 + (tCk-to)' 



where Ck = |/k|. In sublattice A for layer 1, the spin- up 



wave functions are 



1 - A/£;, 



1/2 



(8) 
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The self-consistent equation for the gap is 

4 

n l\ 
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where /i' = /i — Ux/2, is the Heaviside step-function, 
and V-Qz is the volume of the Brillouin zone. The total 
number of electrons (per site) n is related to /i according 
to 



S — 1 



(10) 



At half-filling, n = 1, a; = 0, and ^' = 0. The lower 
two bands are filled while the upper two are empty. Upon 
electron doping, cc > (hole doping, x < 0), ^' abruptly 
changes to the new value > A (/i' < —A). Substitut- 
ing the wave functions wj^^j^ into Eq. (|9|), one obtains 



1 = ^ JdCpoiC) 



i-e(^|M'IA-vW^C + clJ7) 
^/s^ + iC + Cof 



l~&(\fi'\/t-JS^ + {(-(,) 



+ (C - Co)' 



(11) 



where 5 = A/t, = to/t, and po(C) is the dimension- 
less density of states po{C) = J dkS{C — Ck)/VBZ. Equa- 
tion (fTO|) implies 



0(l/^'IA-\/<52 + (C + Co)') + 



(12) 



Solving Eqs. (|lip and ([T^ we obtain A{x) and fJ,{x). 
This can be done analytically if Ao ^ t, to (Aq is the gap 
at zero doping). If Aq is small, the value of ~ Aq is 
also small, and we can omit 0-functions in the first terms 
in Eqs. PT|) and p^ . From these, we derive 

2po(Co)ln(Ao/A) - 2po(Co) asinh(<5C/'5) , 

\x\ - PoiCo)SC, (13) 

where SC = yjiix'f - A^ /t. Solving Eqs. we obtain 



A = AoVM2^, 
/i = Aq [sgn {x) — x/2xc] + Ux/2, 



(14) 
(15) 



where the critical doping Xc ~ Aoto/n^/St^ (the analyt- 
ical expression for Aq in the limit Aq ^ t, to was found 
in Ref. |l3|) . We see from Eq. ((Ti)) that the value of the 
gap decreases with doping, and A = 0, if \x\ > Xc- The 
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FIG. 2: (Color online) (x, U) phase diagram. Solid (red) 
lines show the boundary of the uniform AFM state Xc- For 
large U , our mean-field calculations are not quantitatively 
valid. To emphasize this, the dotted lines plot Xc for U > 
9 eV. The dashed (blue) lines show the boundary of the phase- 
separated state. The inset shows the magnified phase diagram 
for 5 eV < [/ <6 eV. 
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FIG. 3: (Color online) Chemical potential ^ versus doping x 
for the homogeneous state, U — 7 eV [solid (blue) line]. The 
horizontal (red) line shows the Maxwell construction, shaded 
areas are equal: 5*1 — 82- 



curves A(a;) are symmetric for electron {x > 0) and hole 
(x < 0) doping. Next-nearest-neighbor hopping breaks 
this symmetry. However, for the parameters character- 
istic of graphene systems, the asymmetry of A(x) does 
not exceed 1-2%. The critical doping function of 

U is shown in Fig. [51 Strictly speaking, Eqs. p4ll5p are 
not valid for Aq ^ t,io- However, numerical calculations 
demonstrate that Eq. (|14p holds true with very high ac- 
curacy for any ratio of Aq/I. 

In addition to the usual AFM order parameter, more 
exotic possibilities are considered in the literature. For 
example, do ping suppresses the AFM gap, inducing a 
canted state [17| , in which the angle between the magne- 
tization vectors in different magnetic sublattices differs 
from 180°. However, our direct numerical calculations of 
the free energy show that such canted state is unstable 
for any doping. Furthermore, the doped AA-BLG is a 
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typical system with imperfect nesting and, therefore, a 
helical AFM state can be induced in it [l^l- This possi- 
bility will be analyzed below. 

Phase separation and metal-insulator transition. — 
The chemical potential /i versus doping obtained from 
Eqs. (HHIIl]) is shown in Fig.|3l Note that ^^l/^x < 0, if 
|x| is small (this result does not depend on the sign of x). 
In particular, from Eq. ()15p it follows that dfi/dx < 0, if 
U/t < TT^/St/to, which is valid for our choice of param- 
eters. Thus, the compressibility k ~ dx/d^ is negative, 
indicating the instability of the homogeneous phase to- 
ward phase separation upon doping. From Fig. [3l there 
are two stable phases with different doping: = and 
xi > 0. The value of xi can be found using the Maxwell 
construction [l^ , according to which the shaded areas in 
Fig.[3]are equal: = 5*2. The calculated values of xi are 
shown by the (blue) dashed lines in Fig.[2]for different U s. 
For the case shown in Fig. [3l xi < Xc, and the uniform 
system separates into AFM insulator and AFM metal. 
For smaller J7, the situation changes: xi > Xc, and the 
co-existing phases are AFM insulator and paramagnetic 
(PM) metal (see the inset in Fig. [2]). 

If the doped system were to remain uniform, even small 
doping would cause a transition from the insulating mag- 
netic phase to a metallic phase, magnetic or not. How- 
ever, the instability of the uniform phase and the ensuing 
phase separation delays the transition to the conduct- 
ing phase until a finite critical concentration of dopants 
is reached. Because of this phase separation, the doped 
charge segregates into clusters inside the insulating AFM 
matrix. The precise structure of such phase depends on 
a variety of factors: impurities and defects in the sam- 
ple or the substrate, the long-range Coulomb repulsion 
that arises due to local charge- neutrality breaking [20| . 



surface tension at the phase boundaries |20|, |21[, and 
electron-phonon interactions [22|. Charge conservation 
implies that the concentration p of the metallic phase is 
p = |x|/xi. The percolative M-I transition occurs if p ex- 
ceeds some threshold value p* , which is usually about 0.5 
for 2D systems, and the corresponding threshold value of 
doping can be estimated as 0.5 xi. 

Discussion. — The most direct and controllable way to 
switch AA-BLG from AFM insulator to metal is doping 
the system with electron or holes, which could be at- 
tained by using appropriate dopants (e.g., NO2 Ca, 
K I23I ) ■ choosing the substrate and applying a gate volt- 
age |24l.[25| or combining these factors. Our analysis pre- 
dicts that, for interaction and hopping parameters values 
typical for graphene systems, phase separation exists in 
the doping range < x < xi, where xi ~ 0.5—1.5%. 
Thus, the M-I transition occurs at x* 0.25-0.75%. 
For graphene systems, the doping levels ~ 1% are within 
the reach of current experimental techniques such as the 
adsorption of NO2 gas molecules [1, [2^. Moreover, even 
higher dopings, necessary to reach the van Hove singu- 
larity, were achieved [2^. These results suggest that the 



M-I transition we discuss in this paper can be realized 
experimentally. 

As mentioned above, we did not include the helical 
AFM state into our considerations. Such simplifica- 
tion may be justified. Indeed, the helical AFM phase 
is mathematically equivalent to the Fulde-Ferrel-Larkin- 
Ovchinnikov (FFLO) state in superconductors (27l. l28j. 
which is very sensitive to disorder (29j and experimen- 
tally difficult to observe. Further, even if the helical 
state survives disorder, the phase separation and the M- 
I transition remain nonetheless: in such a situation the 
electrons segregate into insulating commensurate AFM 
and metallic helical phases [2^ |3o| with the critical con- 
centration X* being slightly different from the values es- 
timated above. At the same time, the mathematical de- 
scription [l^l of the helical AFM is fairly involved and 
cumbersome. Thus, we believe that at the present stage 
of this research our simplification of the M-I transition is 
warranted. 

The above calculations are restricted to the mean field 
approximation. To what extent the mean field theory 
offers a reliable description of the system? This ques- 
tion was discussed in Ref. |3l| for the usual BCS model 
and for the BCSdike models with finite spin polariza- 
tion in Refs. 
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-3J. It is generally agreed that for weak 
interaction the mean field calculations are accurate in 
these situations. In the intermediate-coupling regime 
the mean field results remain qualitatively correct. Since 
the superconducting systems investigated in these papers 
are mathematically equivalent to the AFM, both doped 
and undoped, we may conclude that our results are at 
least qualitatively correct even for moderately high U. 
Currently, numerical many -body approaches (functional 

and Monte Carlo 



renormalization group 



demonstrated their usefulness for studies of monolayer 
and bilayer graphene. These methods may be used as 
alternatives to the mean field approach. 

If we want to generalize the formalism for finite T, 
we must remember that in 2D at T > no long-range 
AFM order exists. However, the short-range AFM or- 
der survives up to temperatures T*(x) ~ A(x). Indeed, 

we 



following the approach described in Refs. 13L 38-4ol| 
obtain the estimate of the crossover temperature in our 
model r*(x) ~ rMF(a;) ~ 0.6A(x), where rMF(x) is 
the mean- field transition temperature [i^. Thus, the 
crossover temperature is higher than 100 K even if U is 
as small as 5 eV. 

The phase separation can also be destroyed if the 
temperature exceeds a certain threshold value Tps. To 
calculate Tps we have to replace 1 — (/i' — E'k) by 
fi-Ek - /i') - /(-Bk - in Eq. ^ and G-functions 
by the Fermi distributions in Eq. p2|). [/(e) is the Fermi 
distribution function]. Then, we derive fi = /^(x,r) as 
a function of doping and temperature. If T > Tpg, the 
function /x(x, T) increases monotonously with x. Our nu- 
merical analysis shows that Tpg > 100 K, if t/ > 5.5 eV. 
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In conclusion, antifcrromagnetic order, a metal- 
insulator transition, and phase separation are predicted 
for the doped AA-stacked graphene bilayer. These ef- 
fects can be observed at temperatures up to 100 K or 
even higher. 
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